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We analyze dynamic properties of vertical graphene-base hot-electron transistors (GB-HETs) 
and consider their operation as detectors of terahertz (THz) radiation using the developed device 
model. The GB-HET model accounts for the tunneling electron injection from the emitter, electron 
propagation across the barrier layers with the partial capture into the GB, and the self-consistent 
oscillations of the electric potential and the hole density in the GB (plasma oscillations), as well as 
the quantum capacitance and the electron transit-time effects. Using the proposed device model, 
we calculate the responsivity of GB-HETs operating as THz detectors as a function of the signal 
frequency, applied bias voltages, and the structural parameters. The inclusion of the plasmonic 
effect leads to the possibility of the HET-GBT operation at the frequencies significantly exceeding 
those limited by the characteristic RG-time. It is found that the responsivity of GB-HETs with a 
sufficiently perfect GB exhibits sharp resonant maxima in the THz range of frequencies associated 
with the excitation of plasma oscillations. The positions of these maxima are controlled by the 
applied bias voltages. The GB-HETs can compete with and even surpass other plasmonic THz 
detectors. 


I. INTRODUCTION 

Recently, vertical hot-electron transistors (HETs) with 
the graphene base (GB) and the bulk emitter and col¬ 
lector separated from the base by the barrier layers - 
the hot-electron graphene-base transistors (GB-HETs) - 
made of Si 02 and AI 2 O 3 were fabricated and studied P- 
|4]. These HETs are fairly promising devices despite their 
modest characteristics at the present. Similar devices can 
be based of GL heterostructures with the hBN, WS 2 , and 
other barrier layers [SHE]- The history of different ver¬ 
sions of HETs, including those with the thin metal base 
and the quantum-well (QW) base, in which the carri¬ 
ers are generated from impurities or induced by the ap¬ 
plied voltages, as well as HETs with resonant-tunneling 
emitter, is rather long (see, for example, Refs. [9HI6]). 
Figure 1 shows schematically the GB-HET structure and 
its band diagrams. Depending on the GB doping, the 
base-emitter and collector-base voltages, Vb and Vc, and 
the thicknesses of barrier layers separating the base from 
the emitter and collector, We and Wq, respectively, the 
GB can be filled either with electrons or holes [compare 
Figs. 1(b) and 1(G)]. 

In GB-HETs (as in HETs) , a significant fraction of 
the electrons injected from the top n-type emitter con¬ 
tact due to tunneling through the barrier top crosses the 
GB and reaches the collector, while the fraction of the 
electrons captured into the GB can be rather small. The 
tunneling electrons create the emitter-collector current 
Jc- The variations of the base-emitter voltage Vb result 


in the variation of Jc- 

In the case of the GB with a two-dimensional elec¬ 
tron gas (2DEG), the GB-HET transistors can also be 
referred to as the N-n-N heterostructure HETs with the 
GB of n-type. In another case, when the GB comprises 
a two-dimensional hole (2DHG), the GB-HET transis¬ 
tors can also be called the N-p-N heterostructure bipolar 
transistors (HBTs) with the GB of p-type. 

The electrons captured into the GB support the 
emitter-base current Jb- The transistor gain in the 
common-emitter configuration is equal to go = Jc / Jb = 
(1 — p)/p, where p is the average capture probability of 
the capture of electrons into the GB during their transit. 
One of the main potential advantages of GB-HETs is the 
high-speed operation associated with the combination of 
a short transit time (due to the vertical structure), a 
high gain go (because of a low probability, p, of the hot 
electron capture into the GB), and a low GB resistance 
(owing to a high mobility of the carriers in graphene). 
As shown recently p in the graphene heterostruc¬ 
tures with proper inter-graphene barrier layers can be 
very small. Very small values of p in GB-HETs can lead 
to a fairly high gain. The use of the undoped base with 
the induced carriers, leads to the exclusion of the scat¬ 
tering of the hot electrons crossing the GB with the im¬ 
purities and to an increase of the mobility of the carriers 
localized in the GB. 

By analogy with HETs with the QW-base made of the 
standard materials [HI [19], GB-HETs can exhibit res¬ 
onant response to the incoming signals associated with 
the excitation of the plasma oscillations. The resonant 
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an antenna. The net rectified current 6Jc,uj can serve as 
the output current in the GB-HETs radiation detectors. 
The detector responsivity oc 5Due to 
the possibility of the plasmonic resonances, the rectified 
component (the detector output signal) can be resonantly 
large, similar to that in the HET detector [21] and other 
plasmonic THz detectors using different transistor struc¬ 
tures, including those incorporating graphene |22H34j . 

In this paper, develop the GB-HET device model and 
evaluate the GB-HET characteristics as a radiation de¬ 
tector of radiation, in particular, in the THz range of 
signal frequencies. 


II. DEVICE MODEL AND RELATED 
EQUATIONS 


FIG. 1: (a) Structure of a GB-HET with undoped GB and 

GBT -HET band diagrams (potential profile in the direction 
perpendicular to the GB plane) at different relations between 
Vb and Vc'- (b) with GB filled with holes and (c) with GB 
filled with electrons. Arrows show propagation of electrons 
across the barrier layers and capture of some portion of elec¬ 
trons into GB. 


plasma frequencies are determined by the characteristic 
plasma wave velocity s (which increases with the carrier 
density Eg) and the lateral sizes of the GB 2L, while the 
quality factor of the plasma resonances is mainly limited 
by the carrier momentum relaxation time r associated 
with the scattering on impurities, various imperfections, 
and phonons. Due to the specific features of the carrier 
statistics and dynamics in the graphene layers HO], the 
plasma velocity s > vw, where vw — 10® cm/s is the 
characteristic velocity of the Dirac energy spectrum, and 
s can markedly exceed that in the QW heterostructures. 
This promotes the realization of with the plasmonic res¬ 
onances in the terahertz (THz) range even in the GB- 
HETs with fairly large lateral sizes. The possibility of 
achieving the elevated carrier mobilities in the GB can 
enable sharp plasmonic resonances in GB-HETs at room 
temperatures or above. 

The application of the ac voltages ±i5K;/2 at the sig¬ 
nal frequency w, associated with the incoming radiation 
(received by an antenna) to the side contacts to the GB, 
results in the variation of the local ac potential differ¬ 
ence 5ipc^ = 5ipui{x) (the axis x is directed in the GB 
plane) between the GB and the emitter contact. The 
emitter, base, and collector current densities Je, js, and 
jc include the dc components jIe.Oi Jb, 0 j and jc,o (de¬ 
termined by the applied bias voltages Vb and Vc), the 
ac components SjB,un and 5jc,uj (proportional to 

(5V/j, and (due to the nonlinear dependence of the tun¬ 
neling injection current on the local potential difference 
between the GB and the emitter) the rectified dc com¬ 
ponents SjE,uj, and Sjc, (proportional to |<5V/jp, 

i.e., to the intensity of the incoming radiation received by 


We consider a GB-HET with the highly conducting 
tunneling emitter and collector of the n-type and the un¬ 
doped GB (as shown in Fig. I) with a 2DHG induced by 
the applied bias voltages creating the N-p-N structure 
with the band diagram shown in Fig. 1(b). The GB- 
HET device model accounts for the tunneling injection of 
hot-electrons from the emitter to the barrier layer(above 
the barrier top), their propagation across the barrier lay¬ 
ers, partial capture of hot-electrons into the GB, and the 
excitation of the self-consistent oscillations of the elec¬ 
tric potential and the hole density in the GB (plasma 
oscillations). The quantum capacitance and the electron 
transit-time effects are also taken into account. 

We assume that the bias voltages Vb and Vc are ap¬ 
plied between the GL-base and the emitter and between 
the collector and the emitter contacts, respectively. In 
the framework of the gradual channel approximation [35] , 
which is valid if We, Wq “C 2L, the density of the two- 
dimensional hole gas (2DHG)S = S(a:,t) in the GL-base 
and its local potential tp = (p{x,t) (counted from the po¬ 
tential of the emitter) are related to each other as 


_ / p - [ije-Vu p - M/e - Vhi - Vc 

Vtt e \ We Wc 

where k is the dielectric constant of these layers, e is 
the hole charge, Vu is the built-in voltage between the 
contact material and an undoped GL, and p is the 
2DHG Fer mi e nergy in the GB. In a degenerate 2DHG, 
p = hvwV'^^, where h is the Planck constant. The de¬ 
pendence of the right side of Eq. (I) on the hole Fermi 
energy is interpreted as the effect of quantum capaci¬ 
tance [Ml Itt] . 

The ac voltages ±(5V/j/2 and w are applied between the 
side GB contacts connected with an antenna, so that the 
ac potential of the GB 6pui = Spui{x) obeys the following 
(asymmetric) conditions: 



5pu)\x=±L = ±( 514 ,/ 2 . 


( 2 ) 


































3 


The side contacts can serve as the slot wave guide trans¬ 
forming the incoming THz radiation signals being re¬ 
ceived by an antenna into the ac voltage [see Fig. 1(a)]. 
The ac component of the collector-emitter voltage 5Vc,u} 
can also arise due to the ac potential drop across the load 
resistance. However, in the GB-HET with the wiring 
under consideration, 5Vc^ui can be disregarded provid¬ 
ing that the GB-HET structure is symmetrical (see be¬ 
low). For GB-HETs with the degenerate 2DHG, the 
relation between the variations of the hole density and 
the potential, and 5tp^ (the ac components) are, as 
follows from Eq. (1), can be expressed via the net ca¬ 
pacitance per unit area C — CgCquant/quant F Cg')^ 
which accounts for the geometrical capacitance Cg = 
{Ce + Cc) = (K/47r)(lF^^ -I- Wq^) (with Ce oc 
and Cc oc being the geometrical emitter and 

collector capacitances, respectively) and the quantum 
capacitance [3SJ [37] Cq^ant = = 

(2e^/io/7rli^w^), where Eq and p,o are the pertinent dc 
values of the hole density and the Fermi energy. 


e6T,iu = C (3) 


electric field in the emitter barrier, and the emitter dif¬ 
ferential conductance, respectively. The dc hole Fermi 
energy in the GB fiQ obeys the following equation: 
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Wc 


(7) 

The ac electric field components perpendicular to the 
GB plane in the emitter barrier, as well as in the collector 
barrier, are respectively given by 


5Fe,uj 


C 


quant 


Stfa 


{Cquant + Cg) We 


( 8 ) 


SFc^ui = — 


a 


quant 




{Cquant + Cg) Wc 


(9) 


Using the Shockley-Ramo theorem, one can find the ac 
component the electron current densities coming to the 
base and the collector: 


Considering the electron tunneling from the emitter 
to the states above the top of the barrier (through the 
triangular barrier), the emitter electron tunneling current 
density Je can be presented as: 


SjEB,ui = crEB,ujSFE,uj, SjEC, lj = EC,u)5Fe,uj- (10) 
Here 


JE = fE exp 



(4) 


Here F = feh), is the characteristic tunnel¬ 

ing field, A is the activation energy for electrons in the 
contact, m is the effective electron mass in the bar¬ 
rier material, a ^ 1 is a numerical coefficient, Fe = 
{(f — Vbi — file)IWE is the electric field in the emitter 
barrier, and is the maximum current density which 
can be provided the emitter contact, As follows from 
Eqs. (2) and (3), the ac component and the rectified com¬ 
ponent of the emitter tunneling current 6jE,ui and 5 
(for F ^ Fq) are, respectively, given by 


5jE,uj = jE,0 j-,2 ^Pe,uj = 0'eSFe,ui- (5) 
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where v is the drift velocity of the hot electrons cross¬ 
ing the barriers above their tops (which is assumed to 
be constant) and te = WeJv and tc = Wcjv are 
the electron transit time across the emitter and collec¬ 
tor barrier layers. The axis z is directed perpendicu¬ 
lar to the GB plane. Since under the boundary condi¬ 
tions (2), 5ipiu{x) = —Sfuii—x) (see below) and, hence, 
5Fe,ui{x) = —6Fe,ui[x), the net ac currents 


Here jE,o = Jb exp(-F/FE,o)), Fe,o = {Ve - Vu - 
PoIcJIWe-, and (Je = jo{F/F^ q) are the emitter dc cur¬ 
rent density (in the absence of the ac signals), the dc 


^Jeb,uj 


J dx6jEB,ui = CPeB,u 


dxSFE 1 


= 0 , 

(13) 
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/ L pL 

dx6jEC,ui = crEc,u / dx6FE,uj = 0 . 
L J~L 


(14) 

Simultaneously for the rectified components of the dc 
current densities from the emitter to the GB and from 
the emitter to the collector one obtains 


djEB,uj — 


PUeF 


6Fej 


F, 


E,0 


= pT 


(15) 


- {^-p)(XeF 

OjEC,uj — — - . 


SFe^ui 

Fe,o 


(1-p)r|5v5^|2, (16) 


where 


^quant 

{C^^ant + CgY ■ 


(17) 


Consequently, the net rectified components of the dc 
emittter-base and emitter-collector currents are given by 


5Jeb,uj 



dxSjEBM 


pT f dx\S(p^\'^, 


(18) 


L _ i-L 

dxSj EC, UJ ^ i^-p)^ / dx\5ip^\'^, (19) 

L J-L 

respectively. 

The ac hole current along the GB is given by 



, ( 20 ) 

x—L 

where the lateral ac conductivity of the GB aBB,ui is 
given by (see, for example. Refs. 


5J, 


BB,ui — —O'BB, 


dS(Pu 
’ dx 


ie'^po 

t/ . .'/■ 

7Th[LU + t/r) 


( 21 ) 


III. PLASMA OSCILLATIONS AND RECTIFIED 
CURRENT 

To calculate the rectified current components using 
Eqs. (15) and (16), one needs to find the spatial dis¬ 
tributions of the ac potentials in the GB Sip^^ = 5tp,^{x). 

For this purpose we use the hydrodynamic equations 
for the hole transport along the GB [43l |44] (see also 
Refs. [211122]) coupled with the Poisson equation solved 
using the gradual channel approximation [i.e., using 
Eq. (1)]. Linearizing the hydrodynamic equations and 


Eq. (1) and taking into account that the ac component 
of the hole Fermi energy is expressed via the variation of 
their density and the potential, we arrive at the following 
equation for the ac component of the GB potential 
(compare with Refs. [111 HU HI] which should be solved 
with boundary conditions given by Eq. (2): 


d'^Sip^ 

dx'^ 


(oj + iiy)(uj + iv) 


= 0 . 


( 22 ) 


Here s is the characteristic velocity of the plasma waves 
in the gated graphene layers, which is given by s = 
^oc [18], m = oc -s/Eo being 

the so-called ’’fictituous” effective hole (electron) mass in 
graphene layers [2Z], v = {(Xb,uj iWsCg) and v = 1 /t + v, 
where r is the hole momentum relaxation time in the 
2DHG and v = Vyisc + f'rad is associated with the con¬ 
tribution of the 2DHG viscosity to the damping (see for 
example, Ref. m) and with the radiation damping of 
the plasma oscillations. The latter mechanism is associ¬ 
ated with the recoil that the holes in the GB feel emitting 
radiation (the pertinent term in the force acting on the 
holes is referred to as the Abraham-Lorentz force or the 
radiation reaction force [46l [47]). Taking into account 
that the viscosity damping rate is proportional to the 
second spatial derivative, for the gated plasmons with 
the acoustic-like spectrum Vyisc = {( is the 2DHG 

viscosity),and considering that Dyad oc {2e^ [45l - 
ETj), we put V = Xjr F pcP’)-, where p is the pertinent 
damping parameter, and c is the speed of light. If the 
viscosity damping surpasses the radiation damping, one 
can set ?7 ~ 

The characteristic plasma-wave velocity s is deter¬ 
mined by So (as well as the thicknesses of the emitter 
and barrier layers We and Wc) [ 20 ] and, hence, can be 
changed by the variations of the bias voltages Ve and 
Vb- Due to this the characteristic plasmonic frequency 
H = TV sjL can be effectively controlled by these voltages. 

Equation (22) with Eq. (2) yield 


^ s.in{^J{uj + w){uj + iv))x/s] 5W 

sin)-^(w -I- iv){uj + iV)L/s\ 2 

One can see from Eq. (23) that the spatial dependence of 
dpyj is rather complex. In particular, when w approaches 
to nO, where n = 1,2,3,..., this distribution can be os¬ 
cillatory with fairly high amplitude of the spatial oscilla¬ 
tions when n ^ ly. 

Substituting 6piy given by Eq. (23) to Eqs. (15) and 
(16), we arrive at the following: 


djEB,ui — pE 

^ sin[^/{uj + iv){uj + iV)x/s] ‘^\5V^ ^^4) 

sin)-^(w -I- iv){uj + iV)L/s\ 4 
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SjEc,Lu ~ (1 -p)r 

sin[^(g; + w){uj + iv)x/s] \5V^ 
sin[-\/ (w + iv){u} + iv)L/s\ 4 


sjebm — pr 


^dx 

sin[-\/(w -|- + w)x/s] 


J~L 

sin)^/(w -I- w){uj + iv)L/s\ 

4 


(26) 


sJec,uj ^ (1 — p)r 


^dx 

sin[-\/(w -|- + ii')x/s] 


J-L 

sin[-\/(w -I- iz/)(w -I- iV)L/s] 

4 


As follows from Eqs. (24) - (27), SJec,uj = [(1 — 
p)/p]SjEB,ui and 6Jec,uj = [(1 - p)/p]SJeb,cj (as for the 
pertinent dc current in the absence of the THz signals), 
so that SJec,uj ^ ^Jeb,u}- 

Using Eqs. (20) and (24), for the ac current between 
the based contacts we obtain 


SJbb,uj = — 


ie^^o cot)-^ {u! + iv){uj + iv)L/s] 


nhs 


uj + iv 5Vu) 
oj + iv 2 


(28) 


Due to the symmetry of the GB-HET structure and the 
asymmetric spatial distribution of there is no rec¬ 
tified component of the lateral current in the GB-base, 
i.e., 6Jbb,lu = 0. 

At very low signal frequencies w <C v, |i/| ~ AnpaE/K, 
from Eqs. (23) and (22) we obtain 


Sipu: 


swl\y{-/vvx/ s) 
sinh('\/r^L/s) 2 ’ 


(29) 


s Jec, uj — (1 — p)r L 


sinh(2-\/zA/L/s) 
{2\/vPL/ s) 


- 1 . 


sinh^ 


4 


~(i-p)rL 


|<5K; 


(30) 


Since in reality \p\ v, there is an intermediate 
range of frequencies \v\ u: < v. Assuming that 
\v\ <C E, < B, from Eq. (27) we arrive at 



FIG. 2: Spatial distributions of rectified component of the 
emitter-collector current density 5jEC,u) at different frequen¬ 
cies UJ (t = 1 ps, s = 2.5 X 10* cm/s, and L — 1.5 


If the characteristic frequency of the plasma oscilla¬ 
tions D = tts/L ;g> V, the ac GB potential amplitude 
\5(pu:\ can markedly exceed the amplitude of the input 
ac signal JV/, when the frequency is close to one of the 
resonant plasma frequencies nil, where n = 1,2,3,... is 
the plasma resonance index. In this case, the rectified 
emitter-collector current is pronouncedly stratified, i.e., 
its density SjEC,ui is a nearly periodic function of the 
coordinate x. 

Figure 2 shows the spatial distributions of the recti¬ 
fied emitter-collector current density (corresponding to 
the current stratification) calculated for different fre¬ 
quencies using Eq. (25). It is assumed that r = 1 ps, 
s = 2.5 X10® cm/s, and L = 1.5 pm. As seen in Fig. 2, the 
spatial stratification of the current is rather pronounced 
when w is close to the plasma resonant frequencies (the 
frequencies U = 5/6 ~ 0.83 THz and 2U/27r = 5/3 ~ 
1.66 THz): the current exhibits two streams centered 
at \x\/L ~ 0.5 when oj ~ 0.83 THz and four streams 
centered at \x\/L ~ 0.25 when uj ~ 1.66 THz. At the 
frequencies far from the resonances, the spatial current 
distribution becomes weakly nonuniform. 


IV. GB-HET DETECTOR RESPONSIVITY 
A. Current responsivity 


sjEc,u.^{i-p)rL 




2 


(X (1 -p) 



(31) 


The rectified current SJec.ui can be considered as the 
output signal used for the detection of electromagnetic 
radiation (in particular, THz radiation). The current 
detector responsivity using this output signal is defined 
as 
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_ 5Jec,uiH 

~ SI^ 


(32) 


where H is the GB-HET lateral size in the direction along 
the contacts to the GB, is intensity of the incident 
radiation, and S is the antenna aperture. The latter 
is given by S' = A^G/47r [50], where G is the antenna 
gain, = 27r c/lu is the radiation wavelength, and c is 
the speed of light in vacuum. Taking into account that 
1^1 = cS^IStt, where E is the radiation electric field in 
vacuum, and estimating SV^j as 5Vuj = XEjT:, one can 
arrive at 


q\2 T 


327r c 


Tu¬ 


rn 


Gonsidering Eqs. (27), (32), and (33), for the current 
responsivity (in A/W units) we find 


«" = f 


^ sin[-Y/(w + jz/)(w + ih')x/s] 


sin[-y /{w + w)(w + iv)L/s] 


(34) 


where 


8(l-p)rLiJ LH 
' cG ’~^Wl 


(35) 


8(1 -p) jo 


F 

F2 


2 


(36) 


As follows from Eq. (34), at relatively long hole mo¬ 
mentum relaxation times t, the GB-HET responsivity 
exhibits a series of very sharp and high peaks. It is in¬ 
structive that in the GB-HETs under consideration, the 
height of the resonant peaks does not decrease with an 
increasing peak index (as in some other devices using 
the plasmonic resonances). Although some decrease in 
the peaks height with the increasing resonance index at¬ 
tributed to the effect of viscosity and radiative damping 
takes place, the pertinent effect is relatively weak (see 
below). 

In the limiting cases ui u, |l/| and |i7| ^ u, < v, 
corresponding to Eqs. (25) and (26), one obtains 


= I?o — gl? (37) 



respectively. The quantities R and depend on the 
geometrical and quantum capacitances and, therefore, on 
the barrier layers thicknesses: 


Ru (X R (X 
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G2 

^quant 


{Cguant 

1 






I 




■ (39) 


For the doping level of the emitter contact Nu = 
10^® cm“®, assuming that the thermal electron velocity 
vt = 10^ cm/s, we obtain j\. = 1.6 x 10® A/cm^. Setting 
A = 0.2 eV and the effective mass in the barrier layer 
m ~ 2.5 X 10“^® g , we arrive at F ^ 2 x 10® V/cm. 
Setting also Fe^ = 5 x 10® V/cm (to provide the hole 
density in the GB about of E = 10^^ cm“2)^ we find 
jEfi — 2.9 X 10^ A/cm^ and <7e — 0.23 Ohm“^cm“^. 
At these parameters, one obtains also Cquant — 2.6 x 
10® cm“^. Using these data and setting in addition 
p <C 1, WeWc/{We + Wc) = 5 nm, k = 4, and G = 1.5, 
from Eq. (36) we obtain p ~ 2 x 10“^ A/W. If L = 1.5 pm 
and FT — 10 pm, Eqs. (35) - (36) yield i? ~ 30 A/W and 
i?o — 20 A/W. These parameters are also used in the 
estimates below. 

In the case of high quality factor of the plasma res¬ 
onances El/iy, the quantities 5Jeb,ui and SJec,ui and, 
hence, the responsivity as functions of the signal fre¬ 
quency w described by Eqs. (27) and (34) exhibit sharp 
peaks at w ~ nU, attributed to the resonant excitation 
of the plasma oscillations (standing plasma waves). The 
peak width is primarily determined by the frequency v. 
Indeed, for El/v 1 Eq. (34) yields 


maxi?,^ ~ An ~ /2H\ ^ 

Using the above estimate for Rq, the peak values of the 
responsivity at U = 5/3 THz and r = 1 ps is approxi¬ 
mately equal to maxA^^ ~ An — 1.33 x 10® A/W. 

Figure 3 shows the dependence of the normalized GB- 
HET current responsivity calculated using Eq. (34) for 
several sets of parameters. As seen, the positions of the 
responsivity peaks shift toward higher frequencies when 
the plasma frequency increases, i.e., when s increases 
and/or L decreases because ft cc s/L (compare the curves 
”1” and ”2”). A shortening of the momentum relaxation 
time T leads to a smearing of the peaks (compare the 
curve ”3” corresponding to r = 1 ps and the curve ”4” 
corresponding to r = 0.5 ps). Figure 4 shows the low¬ 
ering and broadening of the resonance peaks of the GB- 
HET current responsivity with the decreasing momen¬ 
tum relaxation time r described by Eq. (34). As seen, at 
T < 0.3 — 0.4 ps, the responsivity peaks vanish, while at 
r = 1 ps they are fairly sharp and high.. 


B. Voltage responsivity 


The variation of the dc current component SJec,ui 
cause by the ac signals results in a change of the volt- 
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FIG. 3: Normalized respoinsivity versus signal frequency for 
GB-HETs with different parameters: 

1 - r = 1 ps, s = 5 X 10® cm/s, L = 1.5 /rm; 

2 - r = 1 ps, s = 5 X 10® cm/s, L = 2.0 /rm; 

3 - r = 1 ps, s = 2.5 X 10® cm/s, L = 1.5 /rm; 

4 - r = 0.5 ps, s = 5 X 10® cm/s, L = 1.5 /rm. 

These parameters correspond to fi/27r = 5/3, 5/4, 5/6, and 
5/3 THz, respectively. 



FIG. 4: Normalized respoinsivity as a function of signal 
frequency and momentum relaxation time (s = 2.5xl0® cm/s 
and L = 1.5 /rm, i.e., ll/27r = 5/6 THz). 


age drop, <514; = —SVcfi, across the load resistor in the 
collector circuit (see Fig. 1). Considering that this leads 
to an extra variation of the dc emitter and collector cur- 
rents S Je^ = [(yECg/(Cquant+ Cg){WE + Wc)]6Vc,o and 
5Jec,0 = [(1 - p)(^ECg/(Cquant + Cg)]5Vcfi because of 
its dependence on the collector contact dc potential. The 
latter dependence is essentially associated with the effect 
of quantum capacitance. Taking this into account, for the 
voltage responsivity of the GB-HET under consideration 
RZ = SViuH/SItu = 6VcfiH/SIo) we obtain 



sin[-^(w -I- ii'){uj + iE)x/s] 

L J-L 

sin)-^ {uj + w){oj + iV)L/s] 


(41) 


Here 


= R-r 


re 


1 + 


Ca 


2LH 


{Cquant + Cg) {We + Wq) 


tccte 


s (42) 


where rc is the load resistance. It is instructive that due 
the the absence of the ac collector current (i.e., the ac 
current through the load resistor), associated with the 
GB-HET structure symmetry and the asymmetry of the 
applied ac signal voltage and the ac potential spatial dis¬ 
tribution along the GB) the RG-factor of the voltage re¬ 
sponsivity is independent of the signal frequency. 

For We = Wc = W, Eqs. (38) and (40) yield 


T>V ^ n __ 

“ r aErc{LHlW^)W ' 

{1 + 271 CquantW/ k) _ 

At rc < (1 + 2ECquantWlK){WlLHaE) 
Eq. (43) yields the obvious formulas: 


(43) 


rc. 


Ro = Rorc (44) 

at low frequencies, and 


R 


V 

n 



(45) 


at the plasma resonance a; = H. 

At rc > (1 + 27TCquantW/ k){W/LH aE) = rc, from 
Eq. (43) we obtain, respectively. 


Cquant^^ \ 

^ ~ ° CeW \ K )' 


(46) 


and 


3i?o (2nV{W^/LH) ( 27TCquantW 

For r = 1 ps, s = 5 X 10®cm/s, and L = 1.5 /im (as for the 
curve ”1” in Fig. 3), so that H = 5/3 THz, as well as ce = 
0.23 A/V cm, Cquant = 2.6 X 10® cm“^ (see the estimate 
in the previous subsection), and k = 4 from Eq.(47) we 
obtain the following estimate: Rq ~ 2 x 10®V/W. For 
the above parameters one obtains rc — 150 Ohm. Even 
at smaller rc, the voltage responsivity can be fairly large. 
Setting rc = 5—10 Ohm, we obtain Rq ~ 100—200 V/W 
and Rq ~ (6.6 — 13.2) x 10® V/W, respectively. 
























































Considering Eqs. (42) and (39), we find the following 
dependences of the current and voltage responsivities 
and on the emitter and collector barriers thickness 

UJ 

W (at not too large rc): 


R^ = R^rc oc 


4_2(^2 

^ quant 


k2(1 + 27r CquantWIn)'^ ' 


(48) 


According to Eq. (48), and R^ markedly decrease in 
the range W > n/^TrCquant — 10 nm. 


V. DISCUSSION 

As follows from Eqs. (34) and(40), an increase in u 
with increasing frequency u due to the reinforcement 
of the plasma oscillation damping associated with the 
viscosity and the radiative damping, might lead to the 
gradual lowering of the resonant peaks with their index 
n. However, our estimates show that the contribution of 
these two mechanisms to the net damping is small com¬ 
pared to the damping associated with the hole momen¬ 
tum relaxation (collisional damping). Indeed, disregard¬ 
ing the radiative damping and assuming the 2DHG vis¬ 
cosity to be C = 10 cm^/s (i.e., smaller than in the stan¬ 
dard 2DEG and 2DHG in the GaAs based heterostruc¬ 
tures [m m [SO]) and s = (2.5 — 5) x 10® cm/s, we 
obtain rj ^ — 16) x 10“^^ s. Hence, in the fre¬ 

quency range a;/27r < 2 THz (as in Fig. 2), we find 
f'jjisc = < (6.3 — 25.3) X 10® s“^ ^ 1/r. Therefore, 

the heights of the responsivity peaks in Fig. 2 in the curve 
”3” at w ~ 0.8 and 1.6 THz are virtually equal. How¬ 
ever, the peaks corresponding to higher resonances with 
the frequencies in the range 5-10 THz can be markedly 
lowered and smeared, because in this range Vyisc can be¬ 
come comparable with 1/r. For example, for the same 
values of ? 7 , T = 1 ps, and w = 5-10 THz, we obtain 
VviscT — 0.16 - 0.64. 

Equation (42) describes the saturation of the voltage 
responsivity peak value with increasing load resis¬ 
tance rc- This is associated with the effect of the voltage 
drop across the load on the potential drop between the 
emitter and the base and, hence, the hole Fermi energy 
in the GB, which determine the injection current. Such 
an effect is due to the finite value of the GB quantum 
capacitance (see. Refs. misi) - if Cquant tends to infin¬ 
ity, the emitter-base voltage becomes independent of rc, 
and the saturation of the Rq — rc dependence vanishes. 

The THz detectors using a similar operation princi¬ 
ple and InP double heterojunction bipolar transistors 
(DHBTs) were recently fabricated and studied experi¬ 
mentally [5TVl53j . The estimated responsivities of the 
DHBTs in question for the non-resonant detection regime 
are somewhat smaller but of the same order of magni¬ 
tude than those given by Eq. (43) and the pertinent esti¬ 
mates. However, the experimental values of the respon¬ 
sivity are much smaller than the values predicted above 
for the resonant detection [see Eqs. (44) and (46) and 


the estimates based on these equations]. Apart from the 
parasitic effects and the absence of any spatial coupling 
antenna, this can be attributed to the doping of the base 
in the InP-DHBTs, which inevitably leads to relatively a 
shorter hole momentum relaxation time r compared to 
that in the GB (where the 2DHG is induced by the ap¬ 
plied voltages). Possibly, a higher probability of the hot 
electron capture into the InGaAs base (that had a rela¬ 
tively large thickness of 28 nm) in the DHBTs in com¬ 
parison with the GB-HETs can be an additional factor. 

The GB-HET current and voltage responsivities are 
determined by several characteristics: the characteris¬ 
tics of the tunneling emitter, geometrical characteristics 
of the GB-HET structure, materials of the emitter as 
well as the emitter and collector barrier layers, and ap¬ 
plied bias voltages. The diversity of these factors enables 
the optimization of the GB-HETs operating as resonant 
plasmonic detectors, in particular, an increase in the re¬ 
sponsivity in comparison with the values obtained in the 
above estimates. The resonant plasmonic THz detectors 
can based on not only the GB-HET structure shown in 
Fig. 1(a) (with an extra antenna connected to the GB 
side contacts), but also based on lateral structures with 
the GB contacts forming a periodic array. 

The comparison of the GB-HETs [iHl] and InP- 
DHBTs [511155] with the GB-HETs under consideration 
highlights the following advantages of the latter: (i) a 
longer momentum relaxation time of holes t in the GB; 
(ii) a higher plasma-wave velocity s that enables higher 
resonant plasma frequencies; (iii) a smaller capture prob¬ 
ability of hot electrons into the GB and, consequently, 
larger (or even much larger) fraction of the hot electrons 
reaching the collector; (iv) coupling the incoming THz 
signal to the GB resulting in the absence of the ac cur¬ 
rent in the emitter-collector circuit and prevanting the 
RG effects usually hindering the high-frequency opera¬ 
tion. 


VI. CONCLUSIONS 

We developed an analytical model for vertical het¬ 
erostructure HETs with the GB of the p-type sandwiched 
between the wide-gap emitter and collector layers and the 
N-type contacts. Using this model, we described the GB- 
HET dynamic properties and studied the GB-HET oper¬ 
ation as detectors of THz radiation. The main features 
of the GB-HETs are high hole mobility in the GB, low 
probability the capture of the hot electrons injected from 
the emitter and crossing the GB, and the absence of the 
collector ac current. These features enable pronounced 
voltage-controlled plasmonic response of the GB-HETs 
to the incoming THz radiation, high hot-electron injec¬ 
tion efficiency, and the elimination of the RC-limitations 
leading to elevated the GB-HET current and voltage re¬ 
sponsivities in the THz range of frequencies, particularly 
at the plasmonic resonances at room temperature. This 
might provide the superiority of the GB-HET-based THz 
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detectors over other plasmonic THz detectors based on 
the standard heterostructures . Thus, the THz detectors 
based on the GB-HETs can be interesting for different 
applications. 
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